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This document describes subroutines which have been 
written to facilitate use of the IBM 1627 plotter associated 
with an IBM 1800 computer. These subroutines have been 
implemented under the Multiprogramming Executive Operating 
System (MPX) and all information pertains to the MPX versions. 

The programs in the first section are designed to permit 
graphic presentation of data with a minimum of programming 
effort. Most standard types of plots can be produced by 
calling one of the programs in this section. Options are 
included for linear or logarithmic tickmarks or grids, multi- 
color plotting, the tracing of several curves or sets of 
data points on the same base, and histogram plotting. 

The subroutines of the second section provide the basic 
tools for those who wish to design their own special purpose 
plots. Included are routines for controlling pen movement 
in terms of inches or user defined units, writing alphabetic 
or numeric text, plotting special symbols, changing the scale 
of plots, and isometric presentation of three dimensional data. 

The third section describes special subroutines used by 
the programs of the first two sections. These are included 
primarily for informational purposes since they will be of 
relatively little utility to the average user. 






The fourth section discusses certain other available 
plotting routines. 

Three appendices contain tables giving the core storage 
requirements of each routine, a list of the subroutines 
required by each, and listings of the routines themselves. 

The table of contents has been designed to give the user 
a quick overview of what is available by including a brief 
description of the purpose of each routine. 

Notation 

The standard FORTRAN conventions for variable names are 
used throughout this document. Variables whose names begin 
with the letters I through N are one-word integers, and those 
names begin with the letters A through H or 0 through Z are 
single precision (two-word) floating point numbers. 

In the FORTRAN calling statements for the various sub- 
routines, the parameter lists may contain either variable 
names or appropriate numerical values. The following two 
calling sequences, for example, are equivalent. 

(1) X=10.0 
Y=3 . 32 
1=1 

CALL MOVE(X, Y, I) 

(2) CALL MOVE (10 .0,3,32,1) 



1 . One-call plot routines 


The routines of this section have been designed to con- 
struct the most commonly used types of plots with a minimum 
of programming effort. 

H1ST0 CALL HISTO (Y,N, XSIZ , YSIZ , YL , YU, IBAR) 

Subroutine HISTO plots N values of floating point array 
Y in the form of a histogram inside a rectangular frame of 
XSIZ by YSIZ inches. The N bars are of equal width. YL 
and YU are respectively the user assigned lower and upper 
bounds of the data. Bars for values which exceed these limits 
will only be plotted to th<? appropriate boundary, and a 
message noting the error will be printed on the typewriter 
(cf. PLERR.) 

If YL<0<YU, a horizontal line at Y = 0 will be drawn across 
the plotting area. The bars of the histogram will be drawn as 
if extending between Y = 0 and Y = Y^ whether Y = 0 lies 
within the plotting area or not. (see figure 1 (a).) If 
IBAR = 1 the individual bars will be separated by lines ex- 
tending towards Y = 0 (figure 1 (b)), if IBAR = 0 these lines 
will be omitted (figure 1(c).) 

At the end of plotting the origin is reset at a point even 
with the bottom boundary of the last plot and two inches to 
the right of its right edge. This permits successive calls 
to HISTO to be plotted on separate areas. 




QUICK CALL QUICK (X, Y, N, XS, YS, XL , XU, YL , YU, ISYM , SS, LINE , IR , IP ) 

Subroutine QUICK plots N points whose coordinates, in 
user's units, are stored in floating point arrays X and Y. 

The plot is contained in a rectangular area of XS by YS 
inches. XL and XU are respectively the user -assigned lower 
and upper boundaries of the X data, and YL and YU those 
of the Y data. X or Y values which exceed the assigned 
limits will be plotted on the appropriate boundary, and an 
error message will be printed on the typewriter (see PLERR.) 

If ISYM is a number between one and eight the corresponding 
symbol (see subroutine SYMBL) will be plotted at each point. 

If ISYM is zero no symbol will be drawn. SS is the size of 
the symbol, if any, in inches and must be an integral multiple 
of 0.04” and must not be less than 0.04". If LINE equals 
zero a line will be traced between consecutive points, if 
equal to oue no line will be drawn. 

If IR equals zero, indicating an initial plot, a rectangular 
frame of XS by YS inches will be drawn before plotting the 
data, while a value of one will cause the data alone to be 
plotted over the area. If either the X or Y axes fall within 
the plotting area (i.e., if XL <0 < XU or YL <0 < XU,) they 
will be drawn. For example the X axis (Y = 0) has been 


plotted in figure 2. If tickmarks or full line divisions 











are desired subroutine GRID may be used, followed by one 
or more calls to QUICK with IR = 0. 

A value of one for IP will cause a pause before plotting 
for paper adjustment or pen change. Console START must be 
pressed to continue execution. If IP is zero no pause will 
be made. 

If two or more separate plots are desired subroutine MOVE 
may be used to reposition the pen. The next call to QUICK 
will then reestablish the origin at the new pen location. 

The plot shown in figure 2 was produced by the following 
two calls to QUICK, where XI and Y1 contain the coordinates 
of the 100 points determining the solid line, X2 and Y2 the 
15 points marked with a triangular symbol, and XL1, XU1, etc. 
define the ranges of the corresponding data. 

CALL QUICK (XI ,Y1,100,5.0,3.0,XL1,XU1,YL1,YU1,0,0.0,0,0,0) 
CALL QUICK (X2 ,Y2,15,5.0,3.0,XL2, XU2 , YL2 , YU2 ,5,0.08,1,1,0) 
G RID CALL GRID (XSIZE , YSIZE, IXTYP, IYTYP,XINC, YINC) 

GRID draws a rectangle of XSIZE by YSIZE inches, sub- 
divided by tickmarks or full-line divisions as explained in 
Table 1. IXTYP and XINC specify the type and interval of 


division along the X axis, while IYTYP and YINC specify type 






and interval on the Y axis. 


■*.' * 




IXTYP (IYTYP) XINC (YINC) 


1 Tickmarks inside rectangle 

Interval in inches between 

equally spaced tickmarks 
or lines. 

2 Tickmarks outside rectangle 

3 Full lines across rectangle 

4 Tickmarks inside rectangle 

Length in inches of one com- 
plete logarithmic cycle (1 
to 10) of tickmarks or lines. 

5 Tickmarks outside rectangle 

6 Full lines across rectangle 


Table I. 

All types of divisions can be used on either axis, and the 
two axes are independent. The program listed below produced 
figure 3, (a) to (f) in succession. These are only six of 
many possible combinations. 


// JOB 27 JUL 71 23.506 HRS 
// FOR TEST 27 JUL 71 23.506 HRS 
*ONF WORD INTEGFRS 
*L I ST SOURCE PROGRAM 
*NONPROCFSS PROGRAM 

CALL GR ID(2.0 . 3.0* 1 . 4.0 .2. 2.0 ) 
CALL MOVE (3.0*0.0* — 1 ) 

CALL GRID(2. 0.3.0. 1*2.0. 2*0. 2) 
CALL MOVF ( 3.0 ,0 .0 .-1 ) 

CALL GRID(2. 0.3.0* 1*3*0. 2*0.2) 
CALL MOVE(-6.0.4oO,-l ) 

CALL GRID(2. 0*3. 0*3. 3*0. 2*0.4) 
CALL MOVE ( 3.0 .0 .0.-1 ) 

CALL GRID(2.0»3.0»3.6»0.2.2.0) 
CALL MOVE (3.0. 0.0,— 1) 

CALL GRID! 2.0. 3.0. 6. 6* 2.0 * 1 .0 ) 

CALL EXIT 

END 


2 . 
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1 1 . General Purpose subroutines 

The subroutines described in this section are designed 
to facilitate the creation of special purpose plots. It 
should be noted that the routines SCAL and FACTR affect the 
operation of the other subroutines of this section. 


SCAL CALL SCAL (XL , XU, XSIZE, YL , YU, YSIZE) 


Subroutine SCAL must be called before the routines PLOT 
or PLOTA are used. This subroutine sets the origin at the 
current pen location and defines a plotting area in terras 
of user units. The location of the origin may later be 


changed if desired by use of the subroutine MOVE. Such a 
change will not affect the relation of user units to inches. 

XL - Value of X at origin in user units. 

XU - Values of X at XSIZE inches from the origin in user's 
units . 

XSIZE - Extent of the basic plotting area along the X 
axis in inches. (Note: This definition and 

that of YSIZE do not prevent plotting outside 
of the basic area defined; i.e., no comparison 
of the data with limits is made in the routines 
of this section.) 

XL - Value of Y at origin in user's units. 

YU - Value of Y at YSIZE inches from the origin in user's 
units. 
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YSIZE - Extent of the basic plotting area along 
the Y axis in inches. 

Example Given a set of X data with values ranging 
between 200 and 400, and a corresponding 
set of Y values ranging between -0.5 and 0.5 
to be plotted in an 8” x 10” area, the initia- 
lizing call would be: 

CALL SCALE (200.0,400.0,8.0,-0.5,0.5,10.0) 

FACTR CALL FACTR (XFAC, YFAC) 

After a call to FACTR all pen movements in the X direction 
will be multiplied by the scaling factor XFAC, and all Y 
movements by YFAC. A call to FACTR affects all subroutines 
of this and the preceding section; i « , HTSTO, QUICK, GRID, 
PLOT, PLOTA, MOVE, GRID, SYMBL, ALPHA, DECD, TRAN, CIRCL, 
and PLOT3. Where successive calls are used, XFAC and YFAC 
are always relative to the original, or true, scale, rather 
than being cumulative. A call of CALL FACTR( 1 .0 , 1 . 0) will 
always reestablish the original scale. 

PLOT CALL PLOT(X,Y,I) 

Subroutine PLOT moves the pen along a straight line to 
the point whose coordinates are X and Y in user's units. The 
pen will be down (i.e., a line will be drawn) if I = 0 and will 
be raised if I “ 1. 
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PLOTA CALL PLOTA(X,Y,N) 

X and Y are Fortran arrays of coordinates expressed in 
user's units. Subroutine PLOTA traces a segmented line, 
starting at the current pen location ana moving to each 
of the N points (X(I) , Y(I) , 1=1 ,N) in succession. 

MOVE CALL MOVE(X,Y,I) 

Subroutine MOVE mo es the pen to a point (X,Y) inches 
from the origin. If I = 0 the pen will be lowered (i.e., a 
line drawn), if I = 1 the pen will be raised, and if I = -1 
the pen will be raised and the origin reestablished at 
the new pen location. 

TGRID CALL TGRID (XSIZE, YSIZE, XINC, YINC) 

TGRID draws a rectangle of XSIZE by YSIZE inches, with 
tickraarks at successive intervals of XINC inches from the 
origin along the boundaries parallel to the X axis, and at 
intervals of YINC inches along those parallel to the Y axis. 

A rectangle or frame without tickmarks may be drawn if desired 
by setting XINC equal to XSIZE, and YINC to YSIZE. 

SYMBL CALL SYMBL (I , S) 

Subroutine SYMBL causes a symbol selected by the value 
of I to be plotted, centered at the current pen location. 

S is the diameter of the symbol in inches. S must be an 
integral multiple of 0.04", otherwise loss of origin may 
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occur. The symbols corresponding to various values of I 
are shown in figure 4. If I is less than 1 or greater 
than 8 a return will be made to the calling routine without 
a symbol being plotted, and an error message will be 
printed on the typewriter (see PLERR) . 

I 12345678 

symbol X + aO^t?0 > 

Figure 4. SYMBL output. 

A LPHA CALL ALPHA (N, IA,HGHT, ANGLE) 

Subroutine ALPHA plots a string of alphanumeric or 
special characters with user specified size and angle to 
the X axis. N is the number of characters to be plotted. 

IA is the name of a FORTRAN integer array containing the 
code numbers corresponding to these characters (see table 
II.) HGHT is the character height in inches. The minimum 
value of HGHT is 0.07" and best results are obtained if 
HGHT is an integral multiple of this figure. ANGLE is 
the clockwise angle between the line of characters and 
the plotter X axis expressed in radians. 

Character plotting commences at the current pen position 
when the call to ALPHA is made, and the pen is returned to 
this position at the end of plotting. As an example, the 
coding shown in figure 5 would cause the message "TIME 1" 


to be plotted in 0.21" letters at an angle of 45 deg. 
to the X axis. 


character 

code 

character 

code 

character 

code 

A 

1 

Q 

17 

5 

32 

B 

2 

R 

18 

6 

33 

C 

3 

S 

19 

7 

34 

D 

4 

T 

20 

8 

35 

E 

5 

U 

21 

9 

36 

F 

6 

V 

22 

blank 

37 

G 

7 

W 

23 

J 

48 

H 

8 

X 

24 

• 

40 

I 

9 

Y 

25 

+ 

39 

J 

10 

Z 

26 

- 

45 

K 

11 

0 

27 

* 

41 

L 

12 

1 

28 

/ 

46 

M 

13 

2 

29 

= 

51 

N 

14 

3 

30 

( 

38 

0 

15 

4 

31 

) 

42 

P 

16 






Table II. ALPHA arguments 

If an illegal character code (i.e., less than zero or 
greater than 51,) is found in the array, an error message 
will be printed (see PLERR) and that character ignored. 


DIMENSION IA (6) 

DATA I A/20 , 9, 13 , 5, 37, 28/ 

CALL ALPHA (6 , IA , 0 . 21 , 0 . 7854) 

Figure 5. Call to ALPHA 
DECD CALL DECD (F , HGHT , ANGLE) 

Subroutine DECD causes the integral portion of the 
floating point number F to be plotted. This integral portion 
may be no larger than 6 digits. If it contains less than 
six significant digits, leading zeroes will be added. HGHT 
is the character height in inches and should be an integral 
multiple of 0.07”. ANGLE is the clockwise angle between 
the number and the X axis, expressed in radians. 

TRAN CALL TRAN(N, I A, HGHT, ANGLE; 

Subroutine TRAN causes the EBCDIC contents (i.e., "A" 
format) of the Fortran integer array IA to be plotted. This 
subroutine makes it unnecessary for the programmer to code 
his message numerically as described under subroutine ALPHA. 

N is the number of characters, and must not exceed 64. HGHT 
is the character height in inches and should be an integral 
multiple of 0.07". ANGLE is the clockwise angle between the 
message and the plotter X axis, expressed in radians. 

CIRCL CALL CIRCL (X, Y, R) 

CIRCL draws a circle whose radius in inches is R, 
centered at the point whose coordinates in inches are X and Y. 
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SCAL3 



CALL SCAL3 (XL , XU , XSl Z , YL , YU , YSIZ , ZL , ZU , Z SIZ , XANG , ZANG ) 


Subroutine SCAL3 performs the same function for PLOT3 
which SCAL does for PLOT, and must be called before PLOT3 
is used. It establishes the origin at the current pen 
location, and defines scaling factors which permit the 
isometric representation of a 3-dimensional point to be 
located by giving its three coordinates in user units. 

Figure 6 shows the relationship of the various parameters 
to the coordinate system. 

XL, YL,ZL - Values of X,Y,Z at origin in user units. 

XU , YU , ZU - Values of X,Y,Z at XSIZ, YSIZ.ZSIZ inches 
from origin in user units. 

XSIZ, YSIZ, ZSIZ - Extent of the basic plotting area along 


the X,Y,Z axes in inches. (Note: These 

definitions do not prevent plotting out- 
side of the basic area defined; i.e., 
no comparison of the data with any limits 
is made in these routines.) 

XANG, ZANG - Angles in degrees which the plotted represen- 
tations of the X and Z axes make with the 
plotter X axis (see figure 6.) 

The location of the origin may be changed by a new 
call to SCAL3 after the pen has been moved to the desired 





location. 


. ■ ' ' ■ 





Figure 6. SCAL3 parameters 

Note that this group of 3-dimensional plot routines 

» 

does not share the same origin or coordinate system as the 
other plot routines, and they cannot be used together with- 
out extreme care . 

PLOT 3 CALL PL0T3(X, Y,Z, I) 

Subroutine PL0T3 moves the pen to the isometric represen- 
tation of the 3-dimensional point (X,Y,Z) in the coordinate 
system initialized by SCAL3 (see figure 7.) If I = 0 the pen 
will be down (i.e., a line will be drawn) and if I = 1 the 
pen will be raised. 
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Figure 7. Isometric representation of point (X,Y,Z) 

III . System subroutines 

The subroutines described in this section are used 
by those of the previous sections in performing their func- 
tions. They are included mainly for the sake of completeness 
since the majority of programming problems can be solved 
more simply and flexibly with the previously described 
routines. 

PLERR CALL PLERR(N) 

Subroutine PLERR causes an error warning message to 
be printed on the typewriter. The message format is shown 
in figure 8. The integer N is printed in hexadecimal form 
as the first of the group of six numbers (0063 or 99 10 in 
figure 8) and indicates the source of the message. The 



following numbers indicate errors in certain routines: 

1 - Illegal parameter in call to SYMBL . Call is 

ignored . 

2 - Data out of specified limits in HISTO. Data bar 

will be plotted only to the boundary. 

3 - Data out of specified limits in QUICK. Data point 

will be plotted on boundary. 

7 - Illegal character in call to ALPHA. The individual 
character will be ignored. 

00.421 

PLOT ERROR NO. ****** 0063 0000 0000 0000 0000 0000 // 

Figure 8. PLERR message 

PIKUP CALL PIKUP (XF, YF) 

A call to PIKUP sets XF and YF equal respectively to 

the current X and Y expansion factors (cf. FACTR.) 

BRING /BITS CALL BRING CALL BITS 

DC addr DC n 

DC bit 1 S n < 16 

The two e r ;ry subroutine BRING/BITS allows the retrieval 

of data which has been packed in contiguous groups of one 

to sixteen bits. Calls may only be made from assembly 

language routines. 

The subroutine is initialized by a call to BRING, where 
'addr' is the absolute word address and 'bit' the bit number 




at which the data retrieval process is to begin. The first 
call to BITS, followed by a positive data constant less than 
or equal to sixteen, will bring that number of bits into the 
accumulator, starting at the specified address. Each 
successive call fetches the next sequential group of bits, 
whose count is determined by the calling parameter. Word 
boundaries are ignored. 

This routine is used by both ALPHA and SYMBL . In both 
of these a pen command is stored as two four bit groups 
which determine the direction and magnitude of the pen move- 
ment, followed by a single bit indicating whether the pen 
is to be up or down. These nine bit command groups are 
packed contiguously in data arrays at the end of each 
routine . 

PLOTX CALL PLOTX 

DC LIST 

LIST DC *-* completion code 

DC /OOOn 

PLOTX transmits a command determined by the value of n, 
the fourth hexadecimal digit of the second control parameter. 
This command may put the plotter on or off-line, raise or 
lower the pen, or cause the pen to be moved one plotter 
increment (0.01") along the X or Y axes or diagonally in a 
combination X and Y movement. Table III lists the commands 
corresponding to each value. For further information on this 
routine see reference 1. 
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u til'll ' * i 


Value 


Plotter command 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 
A 
B 


C-F 


Pen down 
Line segment +Y 
" " +X,+Y 

" " +X 

" " +X,-Y 

t» tf _y 

” " -X, -Y 

" " -X 

" " -X, +Y 

Pen up 

Take plotter off line 
Put plotter on line 
Not used 


Table III. PLOTX commands 

LINE/SKIP CALL LINE (IX, I Y) 

CALL SKIP (IX, IY) 

A call to subroutine LINE causes the plotter to draw 
the closest approximation of a straight line from the current 
pen location to a given point. IX and IY are the integer X 
and Y coordinates of this point in hundredths of an inch, 
with the starting point defined as the origin. A call to 
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SKIP causes the plotter to move to the specified point 
with the pen raised. The straight line approximation 
is determined by the following algorithm. 

Straight line plotting algorithm . 

The plotter moves vertically and horizontally in steps 
of 0.01". There are eight possible plotter movements; vertical 
up or down, horizontal right or left, and the four diagonal 
combinations of horizontal and vertical movement (see 


figure 9.) 




Figure 9. Plotter incremental movements. 


It is assumed that the closest approximation to a / 

directed line segment can be formed using the two plotter 
movements which define the octant in which the segment lies. 

A line segment to be plotted will be defined as a vector 
with two integer components; the number respectively of 
horizontal and vertical 0.01" steps necessary to reach the 
end point of the segment. 

Suppose that P is a line segment with components m and 
n, and that n - m. Then the closest approximation to P 
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is defined by a series of n points. The coordinates of 
the ith point will be i in the direction of the larger 
component, and 

[ 2im + n"l 

2n J 

in the direction of the smaller, where the square brackets 
represent integer truncation of the quotient. The second 
coordinate is the value of the integer nearest to i times 
the quotient m/n. Figure 10 shows the example P = (3,5). 


i 

„ T2i(3) + 5l 

i 1 10 J 

Y i - 1 

1 

[11/10] = 1 

1 

2 

[17/10] - 1 

2 

3 

[23/10] = 2 

3 

4 

[29/10] = 2 

4 

5 

[35/10] = 3 

1 

5 



Figure 10. Approximation of the line P = (3,5) 

Each of the i successive points has its coordinates computed 
in this manner. The pen movement for the ith step is then 
the difference between the coordinates of that point and 
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those of the previous point: 




The successive commands for the example of figure 
would be: 

( 1 , 1 ) ( 0 , 1 ) ( 1 , 1 ) ( 0 , 1 ) ( 1 , 1 ). 

IV. Other Plotting routines 
CCOMP 

The program CCOMP produces a IBM 1637 plot from a 
Cal Comp 77C plotter tape written on the IBM 360/75 J com- 
puter. The plot is the nearest approximation to that which 
would be drawn on the 770 consistent with the difference 
in resolution of the two plotters. Ordinarily very little 
perceptible distortion occurs. 

The calling sequences for this program is given below. 
The fourth card contains the numbers of the first and last 
CalComp files to be plotted, in 216 format. The CalComp 
file numbers are three digits integers, and 999 is generally 
used to end the tape. A detailed description of this program 
will be given in a document now in preparation. 

// JOB 

// XEQ CCOMP 

* CCEND four cards 

1 999 


TGRAF and PGRAF 


TGRAF and PGRAF produce X-Y plots on the 1053 type- 
writer and 1443 printer respectively. Their operation is 
described in detail in reference 2. 
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APPENDIX A. SUBROUTINE STORAGE REQUIREMENTS 


Name 

words 

HISTO 

252 

QUICK 

294 

GRID 

484 

SCAL/FACTR/PI KUP/ 
MOVE/PLOT 

204 

PLOTA 

58 

TGRID 

280 

SYMBL 

180 

ALPHA 

398 

DECD 

122 

TRAN 

192 

CIRCL 

242 

SCAL3/PLOT3 

252 

PLERR 

42 

BRING/BITS 

52 

PLOTX 

277 

LINE/SKIP 

132 

PGRAF 

392 

TGRAF 

390 


APPENDIX B. REQUIRED SUBROUTINES 


Name Subi outines used 

HISTO SUBIN, SCAL,FLD,LDFAC, SUBSC , FADD, PLERR, 

FMP Y , FDI V , SNR , L I NE , FARC , PLOTX 

QUICK SUBIN, SCAL, PAUSE, FLD,LDFAC, SUBSC, SYMBL, 

PLERR, FMP Y, FDI V, SNR, FADD, LINE, FARC, 
BRING, PLOTX 

GRID SUBIN , FLD , SUBSC , FMPY , LDFAC , MOVE , FARC , 

FDI V, SNR, FADD , LINE , PLOTX 


SCAL /FACTR/PI KUF/ 

MOVE/PLOT FMPY, FDI V, SNR, FADD, LINE, FLD, FARC , PLOTX 

PLOTA SUBIN, SUBSC, FMPY, FDI V, SNR, FADD, LINE, FLD, 

FARC, PLOTX 


TGftID 

SYMBL 

ALPHA 

DECD 

TRAN 

CIRCL 


SUBI N , MOVE , FLD , LDF SC , FMP Y , FADD , LI NE , 

FDI V , SNR , FARC , PLOTX 

PI KUP , FLD , FDI V , FMP Y , BRI NG , FADD , L INE , 

PLERR , SNR , FARC , PLOTX 

FSIN , FMPY , FDI V, PI KUP , BRING , SNR , FADD , LINE , 
PLERR , XMDS, FARC , FLD , PLOTX 

SUBIN, FABS, FLD, FAXI , FDI V, FMPY, SUBSC , I STOX, 
LDFAC , ALPHA , FARC , FSIN , PI KUP , BRING , SNR , FADD , 
LINE , PLERR , XMDS, PLOTX 

ALPHA , FSI N , FMPY , FDI V, PI KUP , BRING , SNR , 

FADD, LINE , PLERR, XMDS, FARC , FLD, PLOTX 

SUBIN , FLD , FADD, MOVE , SUBSC , FMPY , FARC , FDI V , 
SNR, LINE, PLOTX 


PLERR 


none 


BRING/BITS none 


PLOTX 


none 


LINE/SKIP 


PLOTX 


t 


Subroutines used 


SCAL3/PLOT3 


FADD, FMPY, DFI V, SNR, FSIN, FLD,PIKUP, 
L I N E , FARC , XMD S , PLOTX 


II 


EBPRT , PRT , IPCNN 
EBPRT , PRT , ITCNN 
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APPENDIX C 


Subroutine Listings 


28 


HW # mm • ***M» ^FIp(M , ' , w* ; WJ*IL l P l ^PW W* * 



* r 1 * iMfc3ISf..4tti 


SUBROUT INF OUICK (X , Y , N , XS I Z . YS I Z , XL . XU , YL , YU , I S YM , SS I Z , L I NE , I R , I P ) 
DIMENSION X ( 1 ) ,Y ( 1 ) 

IF( IP) 16, 16,17 
17 PAUSE 

16 CALL SCAL(XL,XU,XSIZ,YL,YU,YSIZ) 

IF(IR>12,19,12 

19 TT =X S I Z 
IF(TT-5. 0)40,40, 41 

4 1 TT = 5 • 0 

40 CALL MOVE (TT, 0,0,0) 

CALL MOVE (0.0, 0.0,1 ) 

CALL MOVE ( X S I Z , 0.0,0) 

CALL MOVE (XSI 7.VSIZ .0 ) 

CALL MOVE ( 0.0,YSIZ,0) 

CALL MOVE ( 0.0, 0.0,0) 

IF (XL )6, 8,« 

6 IF(XU)8,8,7 

7 CALL PLOT ( 0 .0 , YL , 1 ) 

CALL PLOT ( 0.0 ♦ YU, 0 ) 

8 IF(YL )9, 12, 12 

9 IF (YU) 1?, 1 2, 10 

10 CALL PLOT (XL,0.0,1) 

CALL PLOT ( XIJ, 0.0 ,0 ) 

12 I E =0 

DO 13 I=1.N 
L-LINE 

I F ( 1-1 ) 30’, 30 ,31 

30 L = 1 

31 XP=X ( I ) 

YP=Y ( I ) 

I F ( XP— XU ) 2 1 ,21,20 

20 XP=XU 

GO TO 23 

21 IF(XL-XP >24.24,22 

22 XP=XL 

23 I E= 1 

24 IF ( YP— YU )26,26,25 

25 YP =YU 

GO TO 28 

26 I F ( YL— YP ) 29 , 29 ,27 

27 YP=YL 

28 I F= 1 

29 CALL PLOT ( XP, YP , L) 

I F ( I SYM ) 1 3 , 13,14 

14 call symbl( isym,ssiz ) 

13 CONTINUF 

I F ( IE) 33 , 33,32 

32 CALL PLERR ( 3 ) 

33 CALL MOVE ( 0.0 ,0 .0 , 1 ) 

RFTURN 

END 


9 • 

' 




SUB ROUT IMP GR ID(XLI M,YLI M. T XTYP, I YTYP, X INC ,Y I NC ) 

DTMFNSJ ON FLOG ( R ) , I PEN ( 3 ) , T I CK ( 3 ) 

DATA I PEN/O ,0 , 1/, TICK /O. 1,-0. 1,0.0/ 

DATA FLOG/O . 30 1 ,0 . 477 , 0 . 60 2 , 0 . 699 , 0 . 778 , 0 .8 4 5 , 0 .90 3 , 0 .9S4 ,1.0/ 

C 

C SUBR. EGRID R/18/70 PROGRAMMER - T. CARLFTON 

c this subroutine will draw a rectangle df sizf xlim (in.) by ylim (in.j. 

C this WILL BF DIVIDcD IN THE X AND Y DIRECTIONS AS SPECIFIED BY THE I NTF- 

C GER PARAMETERS IXTYP AND I YTYP RESPECTIVELY ACCORDING TO THE FOLLOWING 

C TABLE. 

C 1 - TICKMARKS INSIDE RECTANGLE DRAWN AT INTERVALS OF XINC(YINC) 

C INCHES STARTING AT ORIGIN, 

C 2 - TICKMARKS OUTSIDE R ECT ANGLE AT SAME INTERVAL AS ABOVE. 

C 3 - FULL LINES DRAWN ACROSS RECTANGLE AT SAME INTERVAL AS ABOVE. 

C 4 - TICKMARKS INSIDE RECTANGLE WITH LOGARITHMIC CYCLE OF XINC(YINC) 

C INCHES STARTING AT ORIGIN. 

C 5 — OUTSIDE TICKMARKS WITH LOGARITHMIC DIVISION AS ABOVE. 

C 6 - FULL LINES WITH LOGARITHMIC DIVISIONS AS ABOVE. 

C 

I 2= ( I XTYP — 1 )/3 
I 1 = 1 XT YP — 3* I 2 
I 4= ( I YTYP-1 )/3 
13 = 1 YTYP — 37 I 4 
T I CK ( 3 ) =YLI M 
XT =T ICK ( t 1 ) 

TICK ( 3 ) =XLI M 
YT =T ICK ( 13) 

P L 1 =0 
L2 = 0 

10 L 1 =L 1+1 
IE( I 2) 14, 1 4, l l 

11 IF(L 1-10 ) 1 3, 12. 12 

12 L 1 = 1 
L2=L2+1 

13 X=XINC*<L2+FL0G;L 1 ) ) 

GO TO 15 

14 X=L1*XINC 

15 IF(X-XLIM) 16, 17, 17 

16 CALL MOVE(X, 0.0,0) 

CALL MOVE ( X , XT , 0 ) 

CALL MOVE (X,0.0,IPEN(I 1 ) ) 

GO TO 10 

17 CALL MOVE ( XLI M,0 .0,0 ) 

L3=0 
L4 = 0 

20 L3=L3+ 1 
IF ( I 4 ) 24 ,24,21 

21 IF(L3— 10)23,22,22 

22 L3= 1 
L4=L4+ 1 

23 Y=Y I NC * ( L4+FL0G ( L3 ) ) 

GO TO 25 

24 Y=L3*YINC 

25 IF(Y-YLIM)26.27,27 

26 CALL MOVE < XLI M.Y ,0 ) 

CALL MOVE (XLI M— Y T , Y , 0 ) 

CALL MOVE (XLIM,Y,IPEN(T3) ) 

GO TO 20 

27 CALL MOVE (XLIM, YLIM, 0 ) 

I F ( I 1-3)30,38,38 

30 L 1 =L 1—1 
1 F ( I 2 ) 35 , 35 ,31 

31 IF(L 1 ) 32. 32,34 

32 L 2=L 2— 1 
IF (L2) 38, 33, 33 

33 LI =9 

34 X=X I NC* ( L2+FL0G ( LI ) ) 

GO TO 37 

35 IF(L1 >38,38.36 

36 X=L 1 *X INC 

37 CALL M OVE (X.YLIM.O) 

CALL MOVEtX .YLIM-XT.O ) 

CALL MOVE (X .YLIM. 0 ) 

GO TO 30 

38 CALL MOVE(O.O.YLIM.O) 

IF(I3— 3)40,48, 48 

40 L3=L3— 1 
IF(I4)45, 45 .41 

41 IF(L3 >42, 42, 44 
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42 L4=L4-1 

I F ( L4 >48,43,43 

43 L3=9 

44 Y=Y I NC* ( L4+FL0G ( L3 ) ) 
GO TO 47 

45 IF(L3 >48,48,46 

46 Y=L3*YINC 

47 CALL MOVE (O.O.Y.O) 
CALL MOVE (YT,Y,o) 
CALL MOVE (0,0, Y,0) 

GO TO 40 

48 CALL MOVE (0,0, 0*0,0) 
RETURN 

END 


31 

30 


1 

2 


3 

4 

5 


7 


8 

9 

10 


15 

6 


1 1 
12 


SUBROUTINE H I STO ( Y , N, XSI Z , YS I Z , YL, YU, I BAR ) 
DIMENSION Y ( 1 ) 

FN=N 


CALL SCAL( 0 . ,FN,XSIZ ,YL.YU. YSIZ ) 
TT=XSIZ 

I F (TT —5,0 ) 30 , 30 ,31 
TT=5.0 

CALL MOVE ( TT, 0.0 ,0 ) 

CALL MOVE (0.0, 0.0,1) 

CALL MOVE ( XSI Z , 0.0,0) 

CALL MOVE (XSIZ»YSIZ»0) 

CALL MOVE ( O.O.YSIZ.O) 

CALL MOVE! 0.0, 0.0,0) 

IF(YL) 1,3,3 
I F ( .'U ) 4 , 4 , 2 
Z=0.0 

CALL PLOT (0.,0.« 1) 

CALL PLOT (FN,0.,0 ) 

GO TO 5 
Z=YL 
GO TO 5 
Z=YU 

CALL PLOT (0«,Z«1) 

IE=0 

XP=0. 

DO 6 1=1 ,N 
YP=Y(I ) 

I F ( YP— YU ) 8 ,8,7 

YP=YU 

IE = 1 

GO TO 10 

IF(YL-YP) 10,10,9 
YP=YL 
I E= 1 

CALL PLOT (XP.YP.O) 

XP = I 

CALL PLOT (XP.YP.O) 

IF(IBAR)6,6,15 
CALL PLOT (XP.Z.O) 

CONTINUE 

CALL MOVE (XSIZ+2.,0.0»— 1 ) 

IF(IE) 12, 12,1 1 
CALL PLERR ( 2 ) 

RETURN 

END 
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SCAL 
LOOP 1 

AD2 


AD4 

PLOT 


PROC2 

proc; 

SAVAO 

c 

J 


ENT 

ENT 


SCAL 

PLOT 

ENT 


PIKUP 

ENT 


FACTR 

ENT 


MOVE 

DC 


0 

BSI 


SAVE 

LDX 

1 

-8 

LD 

2 

0 

MDX 

2 

+ 1 

STO 


AD2 

STO 


AD4 

BSI 


LOAD 

LI BF 


FSUB 

DC 


0 

LI BF 


FSTO 

DC 


TEMP 

BSI 


LOAD 

LI BF 


FMPY 

DC 


F100 

LI BF 


FDI V 

DC 


TEMP 

LI BF 


FSTOX 

DC 


J 

LI RF 


FMPY 

DC 


0 

LI BF 


SNR 

LI RF 


FSTOX 

DC 


J+2 

MDX 

1 

+ 4 

MDX 


LOOP1 

SLA 


1 6 

STO 


J 

STO 


J+ 1 

MDX 


EXIT 

DC 


0 

BSI 


SAVE 

BSI 


LOAD 

LI BF 


FMPY 

DC 


C 

LI BF 


FADD 

DC 


C+2 

LI BF 


FMPY 

DC 


FACT X 

LDX 

LI 

J 

BSI 


PROC1 

BSI 


LOAD 

LI BF 


FMPY 

DC 


C + 4 

LIEF 


FADD 

DC 


C+6 

L I BF 


FMPY 

DC 


facty 

MDX 

1 

+ 1 

BSI 


PROC1 

MDX 

2 

+ 1 

LD 

12 

-1 

BSC 

L 

PROC2. Z 

CALL 


LI NE 

DC 


JX 

DC 


JY 

MDX 


EXIT 

CALL 


SKIP 

DC 


JX 

DC 


JY 

MDX 


EXIT 

DC 


0 

LI BF 


FADD 

DC 


fhalf 

LI BF 


I F I X 

S 

1 

0 

STO 

1 

+ 2 

A 

1 

0 

STO 

1 

0 

BSC 

I 

PROC1 

BSS 

E 

0 

BSS 


2 

BSS 


8 

BSS 


2 


LOOP CONTROL. 

GETT ADDR. OF XL (OR YL ) 

STORE ADDRESS 
STORE ADDRESS 
GET XU (OR YU) 

SUBTRACT XL (OR YL) 

STORE DIFFERENCE. 

GET X (OR Y) DIMENSION 

MULTIPLY BY 100 

DIVIDE BY DIFFERENCE 

STORE INCHES/UNIT (X OR Y) 

MULTIPLY BY XL (OR YL ) 
CHANGE SIGN 

STORE ORIGIN CORRECTION 
INCREMENT & LOOP. 

RESET ORIGIN BY ZEROING X 
& Y PEN POSITION 
ACCUMULATORS 

GET X COORDINATE 
MULTIPLY BY INCHES/UNIT 
ADJUST FOR ORIGIN VALUE 
MULT. BY X EXPANSION FACTOR 

GET Y COORDINATE 
MULTIPLY BY INCHES/UNIT (Y) 
ADJUST FOR ORIGIN 
MULT. BY Y EXP. FACTOR 

GET PEN UP OR DOWN COMMD. 
BRANCH IF NONZERO 

DRAW LINE TO POINT 
MOVE TO PT. WITH PEN UP 


U- 

1 

1 

1 

1 

1 

1 

1 

1 


THIS SECTION ROUNDS OFF 
THE FLOAT. PT. INCREMENT 
COUNT & STORES AS PEN COM. 
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Hf 

i 


JX 

DC 


0 



JY 

DC 


0 



F100 

DEC 


100 .0 



FHALF 

DEC 


0.5 



FACTX 

DEC 


1 .0 



FACTY 

DEC 


1 .0 



TFMP 

DC 


0 



SAVE 

DC 


0 

s- 


r 

ST X 

1 

XR 1+1 

s 


M 

STX 

2 

XR 1 + 3 

s 


f 

STX 

3 

XR 1 + 5 

s 

THIS SECTION SAVES STATUS 

if 

STD 


SAVAO 

s 


§ 

STS 


STAT 

s 

AND LOADS THE ADDRESS OF 

• 

LDX 

12 

SAVE 

s 



MDX 

2 

-2 

s 

THE FIRST CALLING PARAMETER 

s 

STX 

2 

AD9S+1 

s 


AD98 

LDX 

I 2 

0 

s 

INTO INDEX REG. 2 

- 1 

LDX 

13 

103 

s 


I 

BSC 

I 

SAVE 

s- 


EXIT 

STX 

2 

AD99+1 

E- 


XR 1 

LDX 

LI 

0 

E 


: J 

* a 

LDX 

L2 

0 

E 

THIS SECTION RESTORES STA- 

| 

LDX 

L3 

0 

E 

TUS & RETURNS TO CALLING 

1 

LDD 


SAVAQ 

E 

PROGRAM. 

STAT 

LDS 


0 

E 


AD99 

BSC 

L 

0 

E- 


LOAD 

DC 


0 

L- 



LD 

2 

0 

L 



STO 


ADI 

L 

FLOAT. PT. NO. IS LOADED 


LIBF 


FLD 

L 

& INDEX REG. 2 INCREMENTED 

ADI 

DC 


0 

L 



MDX 

2 

+ 1 

L 



BSC 

I 

LOAD 

L- 


P I K UP 

DC 


0 

P- 



BSI 


SAVE 

P 



LDD 


FACTX 

P 

THIS SECTION RETURNS THE 


STD 

12 

0 

P 

X AND Y EXPANSION FACTORS 


LDD 


FACTY 

P 

TO THE CALLING PROGRAM. 


STD 

I 2 

1 

P 



MDX 

2 

+ 2 

P 



MDX 


EXIT 

P- 


MOVE 

DC 


0 




BSI 


SAVE 




LDX 

1 

-2 


LOOP CONTROL 

LOOP2 

BSI 


LOAD 


GET COORDINATE 


LD 

LI 

AD80+2 

X- 



STO 


AD8 1 

X 

MULTIPLY BY X OR Y 


LIBF 


FMPY 

X 

EXPANSION FACTOR 

AD8 1 

DC 


0 

X- 



LIBF 


FMPY 




DC 


FIDO 


MULTIPLY BY 100 


LIBF 


FADD 




DC 


FHALF 


ROUND OFF 


LIBF 


I FI X 




S 

LI 

J+2 


SUBTRACT CURRENT PEN LOC. 


STO 

LI 

JX+2 


STORE DIFF. AS PEN COMMAND 


MDX 

1 

+ 1 




MDX 


LOOP2 


LOOP BACK 


LD 

I 2 

0 


GET PEN UP OR DOWN COMMD 


MDX 

2 

+ 1 




BSC 

L 

PR0C3.Z 


BRANCH IF NONZERO 


CALL 


LI ME 



. 

DC 


JX 


DRAW LINE TO PT . 


DC 


JY 




MDX 


PR OCA 



PROC3 

CALL 


SKIP 




DC 


JX 


MOVE TO PT. WITH PEN UP 


DC 


JY 




BSC 

L 

PR0C4. — 




SLA 


16 


RESET ORIGIN BY ZEROING 


STO 


J 


X AND Y PEN POSITION 


STO 


J+l 


ACCUMULATORS. 


MDX 


EXIT 



PR0C4 

LD 


J 

u- 



A 


JX 

uu 


STO 


J 

u 

THIS SECTION UPDATES THE 


LD 


J+l 

u 

X 6 Y PEN LOC. ACCUMULATORS 


A 


JX+1 

u 



STO 


J+l 

u- 
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THIS SECTION RESETS THE 


MDX 

DC 


EXIT 

0 

F— 

RSI 


SAVE 

F 

LDD 

12 

0 

F 

STD 


FACTX 

F 

LDD 

IP 

1 

F 

STD 


FACTY 

F 

MDX 

2 

+ 2 

F 

MDX 


EXIT 

F- 

DC 

DC 

END 


FACTX 

FACTY 



X E. Y EXPANSION FACTORS 


SUBROUTINE TO R I D ( X L , YL , X I , Y I ) 
CALL P I K UP ( X F * YF ) 

XL IM=XL*XF 
YL IM=YL*YF 
XTNC=XI*XF 
YINC =Y I ♦ YF 
NX = XL I M/X INC 
I F (NX )5 » 5, 1 1 

11 DO 1 1=1, NX 

X = I *X INC 

CALL MOVF ( X ,0 .0 ,0 ) 

CALL MOVE(X ,0. 1 ,0 ) 

1 CALL MOVE (X ,0 ,0 ,0 ) 

5 CALL MOVE(XLIM, 0.0,0) 
NY=YLIM/Y INC 

I F (NY )6 ,6, I 2 

12 DO 2 1=1, NY 
Y=I*YINC 

CALL MOVE (XLT M,Y ,0 ) 

CALL MOVF ( XL I M— 0 • I , Y ,0 ) 

2 CALL MOVF (XLT M,Y ,0 ) 

6 CALL MOVEIXLI M.YLI M,0 ) 

IF (NX )7,7 , 13 

13 DO 3 1=1 , NX 

X= ( NX+ 1 —I ) *X INC 
CALL MOVE (X,YLIM,0) 

CALL MOVF (X,YLIM-0. 1 ,0 ) 

3 CALL MOVFIX , YLI M,0 ) 

7 CALL MOVE (0.0, YLI M,0 ) 

I F | NY )B,H, 1 4 

14 DO 4 1=1 , NY 

Y= (NY+1 -I )*Y INC 
CALL MOVE (0.0, Y,0 ) 

CALL MOVF (O.l.Y.O) 

4 CAuL MOVE (0.0, Y,0 ) 

8 CALL M0VE(0.0 ,0.0 ,0 ) 

RETURN 






LETTFR SIZE 


TEMPORARY 
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- 




L5 

LIBF 


I F I X 



STO 

LI 

Y+4 



MDX 

1 

+ 2 



MDX 


L4 



CALL 


BITS 



DC 


1 



BSC 

L 

SKI P.+- 



CALL 


LI NE 



DC 


X 



DC 


Y 



MDX 


L3 



MDX 


EXIT 


ERROR 

CALL 


PLERR 



DC 


Cl 


EXIT 

LDD 


AO 



MDX 

L 

SYMBL.+2 


RESTO 

LDS 


0 


X R 1 

LDX 

LI 

0 


XR3 

LDX 

L3 

0 



BSC 

I 

SYMBL 


AO 

DEC 


0 


SKIP 

CALL 


SKI P 



DC 


X 



DC 


Y 



MDX 


L3 


Cl 

DC 


1 


ce 

DC 


8 


YF 

DEC 


0 

CARD 1 

XF 

DEC 


0 

CARD2 

FI 

DEC 


0 

CARD3 

F04 

DEC 


.0395 


FHALF 

DEC 


0.5 

CARD 8 

Y 

DC 


0 

CARD9 

T 1 

DC 


0 


X 

DC 


0 

CARD 1 

PI 

DC 


ADD 1 


P2 

DC 


ADD2 


AOD1 

BSS 


0 



DC 


/OOOO 

SOI 


DC 


/3A07 

02 


DC 


/40AE 

S03 


DC 


/0E81 

S04 


DC 


/ 1 9 1 4 

S05 


DC 


/ A 1 7B 

S06 

ADD2 

BSS 


0 



DC 


/A A A 2 

TO 1 


DC 


/6AAD 

T 0 2 


DC 


/5C4B 

T 03 


DC 


/5422 

T04 


DC 


/BO 98 

T05 


DC 


/ AD45 

TO 6 


DC 


/2351 

T 07 


DC 


/O A A C 

TOR 


DC 


/2523 

T 09 


DC 


/19C8 

T 1 0 


DC 


/B jO 2 

T 1 1 


DC 


/2B4D 

T 1 2 


DC 


/66B5 

T 1 3 


DC 


/6AR0 

T 1 4 


DC 


/CO 8 A 

T 1 5 


DC 


/D252 

T 1 6 


DC 


/3598 

T 1 7 


DC 


/60 43 

T 1 8 


DC 


/5B39 

T 1 9 


DC 


/I 64C 

T20 


DC 


/50 1 A 

T21 


DC 


/3958 

T22 


DC 


/4A56 

T2 3 


DC 


/AO 15 

T24 


DC 


/I 'i6C 

T25 


DC 


/671 A 

T26 


DC 


/OOOO 

T2 7 


END 
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ENT 


ALPHA 


AO 

DEC 


0 


C 46 

DC 


51 

TOTAL NO, OF CHARACTERS 

FXIT 

LDD 


AO 


RESTO 

LDS 


0 


XR1 

LDX 

LI 

0 


XR2 

LDX 

L2 

0 


XR3 

LDX 

L3 

0 



BSC 

I 

ALPHA 


ALPHA 

DC 


0 



STD 


AQ 



STS 


RESTO 



STX 

1 

XR 1 + 1 



ST X 

2 

XR2+1 



STX 

3 

XR 3+ 1 



LDX 

13 

103 

TRANSFER VECTOR ADDRESS 


LDX 

I 1 

ALPHA 



LD 

I 1 

0 



STO 

L 

LC 



LD 

1 

1 



STO 


Ll + 1 



LDD 

I 1 

+ 2 



STD 

L 

F 3 

LETTER SIZE 


LDD 

I l 

+ 3 



STD 

L 

F2 

LETTER ANGLE 


CALL 


FSIN 



DC 


F2 



LI RF 


FMPY 



DC 


F3 



LI BF 


FDI V 



DC 


FO 7 



L I RF 


FSTO 



DC 


FI 



CALL 


FCOS 



DC 


F2 



L I BF 


FMPY 



DC 


F3 



LI BF 


FDI V 



DC 


FO 7 



LI BF 


FSTO 



DC 


F2 



MDX 

L 

ALPHA, +4 



CALL 


PI KUP 



DC 


XF 



DC 


YF 



SLA 


16 



STO 


XTOT 



STO 


YTOT 


LI 

LD 

L 

0 



MDX 

L 

Ll+1 ,-l 



CMP 


C46 



MDX 


ERROR 



NOP 

BSC 

L 

ERR OR, + 



S 


Cl 

TEMPORARY 


STO 


T 1 



SLA 


1 



A 


T 1 



SLA 


2 



SRT 


4 



A 

L 

PI 



STO 


0 



SR A 


16 



SLT 


4 



STO 


R 



CALL 


BRING 


0 

DC 


0 


R 

DC 


0 


L2 

CALL 


BITS 



DC 


8 



A 


P2 



STO 


S 



CALL 


BITS 



DC 


4 



STO 


T 



CALL 


BRING 


S 

DC 


0 


T 

DC 


0 


L3 

CALL 


BI TS 
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DC 


4 



BSC 

L 

L2 * +- 



S 


Cl 



BSC 

L 

NEXT * + — 



S 


C7 



STO 


X 



CALL 


BITS 



DC 


4 



S 


C8 



STO 


Y 



LI RF 


FLOAT 



LI BF 


FMPY 



DC 


F2 



LI BF 


FSTO 



DC 


F3 



LD 


Y 



L I BF 


FLOAT 



LI BF 


FMPY 



DC 


FI 



LI BF 


SNR 



L I BF 


FSTO 



DC 


F4 



LDX 

1 

-4 


L4 

LD 


X 



L I BF 


FLOAT 



LI BF 


FMPYX 



DC 


Fl+4 



LI BF 


FADDX 



DC 


F3+4 



LI BF 


FMPYX 



DC 


YF+4 



LI BF 


I FI X 



BSC 

L 

NEGA.Z+ 



LI BF 


FADD 



DC 


FHALF 



MDX 


L5 


NEGA 

LI BF 


FSUB 



DC 


FHALF 


L5 

L I BF 


I FI X 



STO 

LI 

Y+4 



MDX 

1 

+ 2 



MDX 


L4 



LD 


XTOT 



S 


X 



STO 


XTOT 



LD 


YTOT 



S 


Y 



STO 


YTOT 



CALL 


BITS 



DC 


1 



BSC 

L 

SKI P.+- 



CALL 


LINE 



DC 


X 



DC 


Y 



MDX 


L3 


SKIP 

CALL 


SKIP 



DC 


X 



DC 


Y 



MDX 


L3 


X TOT 

DC 


0 


YTOT 

DC 


0 


Cl 

DC 


1 


C7 

DC 


7 


CB 

DC 


8 


T1 

DC 


0 


ERROR 

CALL 


PL ERR 



DC 


C7 


NEXT 

MDX 

L 

LC.-l 



MDX 


LI 



CALL 


SKIP 



DC 


XTOT 



DC 


YTOT 



BSC 

L 

EXIT 


YF 

DEC 


0 

CARD 1 

XF 

DEC 


0 

CARD2 

FI 

DEC 


0 

CARD3 

F2 

DEC 


0 

CA9D4 

F 3 

DEC 


0 

CARDS 

F4 

DEC 


0 

CARD 6 





F 07 

DEC 

0.069 



FHALF 

DEC 

0.5 

CARD 

a 

Y 

DC 

0 

CARD9 

LC 

DC 

0 

CARD 

10 

X 

DC 

0 

CARD 

1 1 

PI 

DC 

ADD1 



P2 

DC 

ADD2 



ADD1 

DC 

/oooo 

A! 



DC 

/430B 

A2 



DC 

/3 1 08 

A3 



DC 

/I 4R1 

A 4 



DC 

/5D1 8 

A5 



DC 

/E1B9 




DC 

/I DD1 

A 7 



DC 

/F323 

A08 



DC 

/6270 

A09 



DC 

/2982 

A10 



DC 

/C930 

A 1 1 



DC 

/6343 

A 1 2 



DC 

/2EB3 

A 1 3 



DC 

/283S 

A 1 4 



DC 

/63E7 

A 1 5 



DC 

/41 84 

A 1 6 



DC 

/3A46 

A 1 7 



DC 

/248D 

A 1 8 



DC 

/4B54 

A 1 9 



DC 

/EF5i 

A20 



DC 

/A53C 

A21 



DC 

/5645 

A22 



DC 

/C564 

A23 



DC 

/7678 

A24 



DC 

/6AC6 

A25 



DC 

/F271 

A26 



DC 

/4751 

A27 



DC 

/I 3E8 

A28 



DC 

/C57B 

A029 



DC 

/2 8 44 

A030 



DC 

/7EC8 

A031 



DC 

/F692 

A032 



DC 

/7927 

A033 



DC 

/7*"D8 

A03 4 



DC 

/7E92 

A035 



DC 

/•» 829 

A036 



DC 

/9279 

A037 



DC 

/2789 

A038 



DC 

/4000 

A 0 39 


ADD2 

DC 

/8DD5 

R! 



DC 

/69B0 

B02 



DC 

/78BA 

B03 



DC 

/4794 

B04 



DC 

/31FB 

R05 



DC 

/8CBE 

B06 



DC 

/1EEF 

B07 



DC 

/58DC 

BOB 



DC 

/65F0 

B09 



DC 

/D77A 

BIO 



DC 

/C7A0 

B 1 1 



DC 

/3927 

B1 2 



DC 

/7B45 

B 1 3 



DC 

/E71B 

R 1 4 



DC 

/99D4 

B15 



DC 

/65F4 

B 1 6 



DC 

/418F 

B 1 7 



DC 

/D469 

B 1 8 



DC 

/B0B6 

B 1 9 



DC 

/6B47 

B20 



DC 

/AO 39 

B21 



DC 

/048C 

B22 



DC 

/7F22 

B23 



DC 

/8AB8 

B24 



DC 

/DA06 

B25 



DC 

/AD51 

B2 6 



DC 

/8680 

B27 



DC 

/5E47 

B28 



DC 

/FI 79 

B29 



DC 

/ 1 84C 

B30 



DC 

/7E84 

B 31 



DC 

/31 48 

B32 



DC 

/7CBE 

B33 



39 






1 





DC 

/ 1 0 AD 

B 1 1 4 

DC 

/0R94 

BUS 

DC 

/5666 

B 1 1 6 

DC 

/F389 

B 1 17 

DC 

/82E7 

B 1 1 8 

DC 

/661D 

B1 19 

DC 

/EDI 7 

B120 

DC 

/9C5fi 

81 21 

DC 

/6751 

B1 22 

DC 

/97C 1 

B 1 23 

DC 

/Cl 87 

B1 24 

DC 

/EA91 

HI 25 

DC 

/9624 

B126 

DC 

/6323 

B 1 27 

DC 

/A819 

B1 28 

DC 

/E525 

B1 29 

DC 

/A1 59 

B 1 30 

DC 

/5603 

B 1 31 

DC 

/EPAB 

B1 32 

DC 

/A66D 

B1 33 

DC 

/0A84 

B! 34 

DC 

/4DE3 

B1 35 

DC 

/OF98 

B136 

DC 

/E407 

B 1 37 

DC 

/3F42 

B 1 38 

DC 

/ 1 8B6 

B1 39 

DC 

/4529 

B 1 40 

DC 

/91 A3 

B1 41 

DC 

/0D41 

B142 

DC 

/EB1 7 

B 1 43 

DC 

/9AE0 

B1 44 

DC 

/86A1 

B 1 45 

DC 

/358B 

B146 

DC 

/BD70 

B1 47 

DC 

/4220 

B1 48 

END 





SUB ROUT I NE PLOTA(XA,YA,N) 
DIMENSION XA(1),YA(1) 

DO 1 1=1, N 

CALL PLOT (XA(I),YA(I),0) 

RETURN 

END 


/ 
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ENT 


TRAN 



000 10 

TRAN 

DC 


0 



00020 


STD 

1. 

TEMP 

SAVE STATUS 

00030 


STS 


SAVF 



000 AO 


STX 

1 

SVXR1+1 



00050 


STX 

? 

SVXP1+3 



00060 


STX 

3 

SVXR 1+5 



00070 


SLT 


32 



00080 


LDX 

I 1 

TRAN 

PICK UP 

A DD ( A DD ) OF N 

00090 


LD 

I 1 

0 

LOAD N 


00 100 


SR A 


1 



00110 


SLA 


1 



on i ?n 


SLA 


9 



001 30 


SR A 


9 



00 1 40 


BSC 

L 

ABORT. + 



001 50 


STO 

L 

WDCNT 

STORE N 

IN WDCNT 

00 1 60 


MDX 

L 

WDCNT.-l 



00 170 


LDX 

I ? 

WDCNT 

STORE N 

IN XR2 

00 160 


MDX 

L 

TRAN. + 1 



00 190 


NOP 





00200 


LD 

I 

TRAN 



UU210 


STO 


ADDZ+ 1 



00 220 

ADD? 

LDX 

LI 

0 



00230 

ADD 

LD 

1 

0 

LOAD FIRST EBCDIC CHARA 

002 40 


LDX 

3 

+ 2 



00250 


RTE 


0 



00260 


STO 


EBCD 



00270 

PFPT 

LD 


FBCD 



00280 


S 


D240 



00290 


BSC 

L 

NE X 1 , - 



00300 


LD 


EBCD 



00310 


S 


D225 



00320 


BSC 

L 

NEX 2 , — Z 



00330 


LD 


EBCD 



00340 


S 


D20B 



00350 


BSC 

L 

NEX3.-Z 



00360 


LD 


EBCD 



00370 


c 


P\1 oo 

k/ 1 ’ t- 



00 38 0 


BSC 

L 

NEXT, -2 



00390 


LD 


EBCD 



00 391 


SLA 


10 



00392 


SR A 


10 



00 39 3 


RTE 


2 



00394 


ERA 


2 



00395 


SLT 


2 



00396 


A 


D37 



00397 


MDX 


NEXT 



00398 

NEX 1 

A 


09 



004C0 

NEX2 

A 


D9 



00410 

NEX 3 

A 


D9 



00420 

NEXT 

NOP 





00430 

ADD? 

STO 

L2 

I A 



00440 


MDX 


RTRN 



00450 


NOP 





00460 

RTRN 

SLA 


16 



00470 


SLT 


8 



00480 


STO 


EBCD 



00490 


MDX 

L 

ADD2+1 .-1 



C0500 


MDX 

3 

-1 



0051 o 


MDX 


REPT 



00520 


MDX 

L 

ADD2+ 1 .+2 



00530 


MDX 

1 

-1 



00540 


MDX 

2 

-2 



0 0 550 


MDX 


ADD 



OOSbO 


LDX 

I 1 

TRAN 



00 570 


MDX 

1 

4. 1 



00 580 


NOP 





00 590 


LD_ 

1 

0 



00 600 


STO 


HGHT 



00610 


MDX 

1 

+ 1 



00620 


NCIP 





00630 


LD 

1 

0 



00640 


STO 


ANGLE 



0 0 650 


LD 


HGHT-1 



00660 


A 


WDC NT 



00670 


STO 


HGHT-1 



G0 6h0 


MDX 

L 

WDC NT, +1 



00590 


LDX 

I 3 

103 

TRANS VECTOR ADDRESS 

00700 


CALL 


alpha 



00710 
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ID tC 



DC 


WDCNT 

00720 


DC 


IA 

00730 

HGHT 

DC 


0 

00740 

ANGLE 

DC 


0 

00750 

ABORT 

MDX 

L 

TRAN, +3 

00760 

SVXR1 

LDX 

LI 

0 

00770 


LDX 

L2 

0 

00780 


LDX 

L3 

0 

00790 

SAVE 

LDS 


0 

00800 


LDD 


TEMP 

00810 


BSC 

I 

TRAN 

00820 

TEMP 

BSS 

E 

2 

00830 

EBCD 

DC 


0 

00840 

D208 

DC 


208 

00850 

D22 5 

DC 


225 

00860 

D240 

DC 


240 

00870 

D192 

DC 


192 

00880 

D9 

DC 


9 

00890 

D37 

DC 


37 

00891 

WDCNT 

DC 


0 

00900 

XR3 

DC 


0 

00910 

PETI N 

DC 


0 

00920 

I A 

BSS 


64 

00930 


END 



0094C 


SUBROUTINE DECD( AA.HGHT, ANGLE) 

DIMENSION I A( 6) , ITBLI 10 ) 

DAT A I TBL/ 27,28,29,30,31,32,33,34,35.36/ 
DATA N/6/ 

A=ABS( AA ) 

D03 J= 1 , N 

B=10.0**J 

FINT=IFIX(A/B) 

C= ( A-F INT*B)/( 10.0**( J-l ) ) 

IC=IFIX(C)+1 

J9=N+1-J 

IA( J9)=ITBL(IC) 

3 CONTINUE 

IF( AA)5,6,6 
I A ( 1 ) =45 
C ONT I NUE 

CALL ALPHA ( N , I A , HGHT , ANGLE ) 

RETURN 
END 



SUBROUTINE CIRCL(X,Y,R) 

DIMENSION SSI IP) 

DATA SS/O .0 .0 .0872, 0.173 6- 0.2588,0.342,0.4226.0.5,0 .5776.0.64 28, 
10.7071, 0.766,0.8192,0.866,0.9063,0.9397,0.9650,0.9848,0 .9962,1.0/ 
CALL MOVE ( X ,Y+R ♦ 1 ) 

DO 1 1=1.18 

J= 1 9— I 

1 CALL MOVE ( X+SS ( I+1)*P,Y+SSI J>*R.O> 

DO 2 1=1,18 

J= 19-1 

2 CALL MOVE ( X+SS ( J)#R,Y-SS( I+1)*R,0) 

DO 3 1=1,18 

J= 19—1 

3 CALL MOVE ( X-SS ( I+1)*R,Y-SSI J)*R,0) 

DO 4 1=1,18 


J= 1 9— I 

4 CALL MOVE (X -SSI 


J)*R,Y+SS( I+1)*R,0> 


RETURN 


END 


ENT 

PLERR DC 
STD 
STS 
LD I 
STO 

A DDR LD L 

STO 

CALL 

DC 

DC 

DC 

DMES 
XIO L 
LDD 

RESTO LDS 

MDX L 
BSC I 
MSG 1 DMES 
MSG 2 DMES 
HEX NO BSS 
AO DEC 
END 


PLERR 

0 

AO 

RESTO 

PLERR 

ADOR+1 

0 

HEXNO 
EACPT 
MSG 1 
MSG2 
HEXNO 
//•E 
46 
AQ 
0 

PLFRR , + 1 
PLFRR 
' RPLOT 'E 
ERROR NO. • E 
6 
0 


44 



ENT 


LINE 



ENT 


SKIP 


LINE 

DC 


0 



BSI 


BEGIN 



LDX 

I 1 

LINE 



BSC 

L 

START, +- 



CALL 


PLOTX 



DC 


L I ST 1 



SLA 


16 



MDX 


START 


SKIP 

DC 


0 



BSI 


BEGIN 



LDX 

I 1 

SKIP 



BSC 

L 

START, Z- 



CALL 


PLOTX 



DC 


LI ST2 



LD 


Cl 


START 

STO 


IND 



STX 

2 

XR2+ 1 



STS 


RESTO 



3LT 


32 



LDX 

2 

0 



LD 

I 1 

0 



BSC 

L 

A 1 » — 



MDX 

2 

+ 8 



SLA 


16 



S 

I 1 

0 

ABSOLUTE VALUE 


CMP 


MAXN 



MDX 


A 1 



MDX 


A 1 



LD 


MAXP 


A 1 

STD 


X 



LD 

I 1 

1 



BSC 

L 

A2 » — 



MDX 

2 

+ 4 



SLA 


16 



S 

I 1 

1 

ABSOLUTE VALUE 


CMP 


MAXN 



MDX 


A2 



MDX 


A2 



LD 


MAXP 


A 2 

STD 


V 



MDX 

1 

+ 2 



STX 

1 

BACK+ 1 



CMP 


X 

Y GREATER THAN X 


MDX 


A3 

IF SO, GO TO A3 


MDX 


A 6 



BSC 

L 

XR1 ,+- 


A6 

SRT 


15 

M2=2Y 


STD 


M2 



LD 


X 



SLA 


1 

N2=2X 


STO 


M2 



SRT 


17 

z 

II 

1- 


STD 


T 1 



LD 


COMM1 



MDX 


A 4 


A3 

SLA 


1 

M2 =2 Y 


5~~ 


N2 



SRT 


17 

H 

i— » 

II 

Z 


3'i D 


T 1 



LD 


X 



SRT 


15 

M2 =2 X 


STD 


M2 



LD 


COMM2 


A 4 

SR A 

2 

0 

SHIFT FOR CORRECT VERT/HORIZ COMMAND 


AMD 


MASK 



STD 


T3 



LD 


COMM3 



SR A 

2 

0 

SHIFT FOR CORRECT DIAGONAL COMMAND 


AMD 


MASK 



STO 


T3+1 



SLA 


1 6 

T2=0 


STO 


T2 



LDX 

I 1 

T1 + 1 

INITIALIZE STEP COUNT 

LOOP 

LDD 


T 1 



AD 


M2 

CALC 2IM+M 


STD 


T 1 



D 


N2 

<2IM+N)/2N 


45 


store RFSULT 



STO 


X 


S 


T2 


A 


ADDR 


STO 


A 5+ 1 

A 5 

LD 

L 

0 


STO 


COMMD 


CALL 


PLOTX 


DC 


LI ST3 


LD 


X 


STO 


T2 


MDX 

1 

-1 


MDX 


LOOP 

XR 1 

LDX 

LI 

0 

X R 2 

LDX 

L2 

0 


LDD 


AQ 

RESTO 

LOS 


0 

RACK 

RSC 

L 

0 

REGIN 

DC 


0 


STX 

1 

XR1 + 1 


STD 


AO 


LD 


I ND 


BSC 

I 

BEGIN 

IND 

DC 


-1 

Cl 

DC 


1 

X 

DC 


0 

Y 

DC 


0 

N2 

DC 


0 

COMM 1 

DC 


/7733 

COMM2 

DC 


/51 51 

COMM3 

DC 


/6842 

MASK 

DC 


/OOOF 

T? 

DC 


0 

ADDR 

DC 


T 3 

M A X N 

DC 


/ROOO 

MAXP 

DC 


/3FFF 


PSS 

E 

0 

AO 

DEC 


0 

T 1 

DEC 


0 

M2 

DEC 


0 

T 3 

DEC 


0 

LI ST 1 

DC 


0 


DC 


0 

LI ST2 

DC 


0 


DC 


o 

LI S T 3 

DC 


0 

COMMD 

DC 


0 


END 


GET COMMAND ADDRESS 


HORIZONTAL COMMANDS 
VERTICAL COMMANDS 
DIAGONAL COMMANDS 


PEN DOWN COMMAND 
PEN UP COMMAND 


46 




ENT 

ENT 

BR I NG DC 
ST X 
LDX 
LD 
STO 
LD 
STO 
LD 
STO 
MDX 
MDX 

BITS DC 
STX 
LDX 

A2 LD 

STO 
MDX 
STX 
BSC 
SLA 
STO 

LOOP LDD 
SLT 
STD 
MDX 
MDX 
MDX 

A 1 LD 

STO 
LD 
STO 

NEXT MDX 
MDX 
LD 

XR1 LDX 

BACK BSC 

Cl 6 DC 

T2 DC 

T 1 DEC 

COUNT DC 

NUMBR DC 
ENn 


BRING 

BITS 

0 

1 XR 1+1 
II BRING 
1 0 
Al + 1 

I 1 0 

Tl + 1 
Cl 6 
COUNT 
1 + 1 
A2 
0 

1 XR 1+1 

II BITS 
1 0 

NUMBR 
1 +1 

1 BACK+1 
L XR1,+- 
16 
T 1 
T 1 
1 

T 1 

L COUNT, -1 
NEXT 

L A1+ 1 , +1 
L 0 

Tl + 1 
Cl 6 
COUNT 

L NUMBR, -1 
LOOP 
T 1 

LI 0 
L 0 
16 
0 
0 
0 
0 
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i 
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LDX 

BSC 

ps dc 

LIRF 

DC 

LIRF 

DC 

LIRF 

DC 

LI RF 
DC 

LTRF 

DC 

LIRF 

DC 

BSC 

SAVE DC 
STX 
STX 
STX 
STD 
STS 
LDX 
MDX 
STX 

AD°R LDX 
LDX 
PSC 

EXIT STX 
XR1 LDX 
LDX 
LDX 
LDD 

STAT LDS 
AD99 BSC 
LOAD DC 
LD 
STO 
LIRF 
ADI DC 
MDX 
RSC 
RSS 

SAVAO RSS 
C BSS 

J BSS 

JX DC 

JY DC 
PLOTS DC 
RSI 
CALL 
DC 
DC 
RSI 
LTRF 
DC 

' I RF 
DC 

LIRF 

DC 

LIRF 

DC 

LTRF 

DC 

LT RF 

DC 

RSI 

LIRF 

DC 

LTRF 

DC 

LIRF 

DC 

RSI 

LTRF 

DC 

LIRF 

DC 

LIRF 


1 6 
I PI 
0 

FSTO 

TEMP 

FMPY 

C+10 

FADDX 

C 

FSTOX 

C 

FLD 
TEMP 
FMPY 
C + fi 
I P3 
0 

1 XR 1+1 

2 XR 1+3 

3 XR1+S 
SAVAO 
STAT 

12 SAVE 
2 -2 

2 AD9P+1 
12 0 

13 103 

I SAVE 
2 AD99+1 
LI 0 
L2 0 
L3 0 

SAVAO 

0 

L 0 
0 

2 0 
ADI 
FLD 
0 

2 +1 
I LOAD 
E 0 
2 

16 

2 

0 

0 

0 

SAVE 
PIKUP 
FACTX 
FACTY 
LOAD 
FSTO 
TEMP2 
FMPY 
C+l 2 
FSTO 
TEMP 
FLD 
TEMP2 
FMPY 
C+l A 
FSTO 
TEMP2 
LOAD 
FMPY 
C+4 
FADD 
TFMP2 
FSTO 
TEMP2 
LOAD 
FSTO 
2 

FMPY 

C+R 

FADD 


MULT. RY Z ORIGIN CORR. 

ADD TO X (Y) ORIGIN CORR. 
NFW X (Y) ORIGIN CORRECTION 


MULT. RY Z UNITS/INCH 

S 

S 

S 

S THIS SECTION SAVES STATUS 
S 

S AND LOADS THE ADDRESS OF 
S 

S THE FIRST CALLING PARAMETER 

5 

S INTO INDFX REG. 2 
S 

s 

F 

E 

E THIS SECTION RESTORES STA- 
E TUS F, RETURNS TO CALLING 
E PROGRAM. 

E 

E 

L 

L 

L FLOAT. PT. NO. IS LOADED 
L 6- INDEX REG. 2 INCREMENTED 
L 
L 

L 


GET X COORDINATE 


/ 


a ' 


GET Y COORDINATE 
MULTIPLY BY INCHES/UNIT 








DC 


TFMP 




LI FF 


FADD 




DC 


C + 2 


ADJUST FDR ORIGIN VALUE 


LI FF 


FMPY 




DC 


FACTX 


MULT. FY X EXPANSION FACTOR 


LDX 

LI 

J 




FST 


PR DC 1 




LT FF 


FLD 




DC 


7 




LI FF 


FMPY 




DC 


C 




LI FF 


FADD 




DC 


TEMP+2 




LT FF 


FADD 




DC 


C + 6 


ADJUST FOR ORIGIN 


LI FF 


FMPY 




DC 


FACTY 


MULT. FY Y EXP. FACTOR 


MDX 

1 

+ 1 




FF I 


PRDC1 




MDX 

2 

+ 1 




LD 

T 2 

-1 


GET PEN UP OR DOWN COMMD. 


FF C 

L 

PROC2 ,7 


RRANCH IF NONZFRO 


CALL 


LTNF 




DC 


JX 


DRAW LIME TO POINT 


DC 


JY 




MDX 


EXIT 



proc2 

CALL 


SKIP 




DC 


JX 


MOVE TO PT. WITH PEN UP 


DC 


JY 




MDX 


EXIT 



proci 

DC 


0 

U- 



LI ff 


FADD 

1 



DC 


FHALF 

1 

THIS SFC.T I ON ROUNDS OFF 


LI FF 


I FI X 

1 



S 

1 

0 

1 

THE FLOAT. PT. INCREMENT 


STO 

1 

+ 2 

1 



A 

1 

0 

1 

COUNT f, STORES AS PEN COM. 


STD 

1 

0 

1 



RSC 

I 

PRDC1 

1- 


Fi no 

DEC 


100.0 



fhalf 

DFC 


0.5 



FACTX 

DEC 


1 .0 



FACTY 

DEC 


1 .0 



PAD 

DEC 


0.0! 7A53 


CONVERSION FACTOR RADIANS 

ANGLE 

DEC 


0 



TEMP 

DFC 


0 



TFMD? 

DEC 


0 



7 

DEC 


0 




END 






